Wrapping of a spherical colloid by a fluid membrane 
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We theoretically study the elastic deformation of a fluid membrane induced by an adhering spher- 
ical colloidal particle within the framework of a Helfrich energy. Based on a full optimization of 
the membrane shape we find a continuous binding and a discontinuous envelopment transition, the 
latter displaying a potentially substantial energy barrier. A small gradient approximation permits 
membrane shape and complex energy to be calculated analyticaUy. While this only leads to a good 
representation of the complex geometry for very small degrees of wrapping, it still gives the correct 
phase boundaries in the regime of low tension. 



INTRODUCTION AND MOTIVATION 



Fluid lipid bilayers are one of the key structural el- 
ements of all living cells. On the one hand they parti- 
tion and thereby organize the complex hierarchy of intra- 
cellular biochemical environments. On the other hand 
they provide controlled transport between neighboring 
compartments as well as the extracellular space. These 
transport mechanisms span quite a large range of par- 
ticle sizes, all the way from sub-nanometer ions (which 
cross the membrane via protein channels) up to micron 
sized objects (like bacteria), which are engulfed in large- 
scale membrane deformations occurring during phagocy- 
tosis 0. Generally, such events require metabolic en- 
ergy and are meticulously controlled by the cell. How- 
ever, there are also cases where they happen passively 
as a consequence of generic physical interactions, for in- 
stance a sufficiently strong adhesion between the particle 
about to be transported and the membrane. An impor- 
tant and well studied example is provided by the route 
along which many animal viruses leave their host cells, 
namely, via the wrapping and subsequent pinch-off of a 
pre-assembled viral nucleoprotein capsid at the plasma 
membrane |2|. Other examples, in which membrane de- 
formations following binding to a small object are cru- 
cial, include certain gene transfection systems, in which 
DNA is complexed by positively charged polymers and 
the resulting condensed globule becomes internalized by 
the cell in an adhesion-driven invagination process ; or 
a host of modern biophysical experimental techniques, in 
which for instance microbeads [i| or AFM tips Q are in 
contact with a membrane and become partially or fully 
wrapped. 

These biological examples are complemented by more 
physically oriented experiments on the adsorption of mi- 
crometer sized beads onto model lipid bilayers. Stud- 
ies have for instance focused on the wrapping of a latex 
bead by a vesicle ^ , or on the interaction between sev- 
eral beads via membrane mediated forces [0. For these 
non-flat substrates, it is unfortunately difficult to extract 
detailed information about the membrane shape close to 
contact, especially for small colloids for which the bend- 
ing contribution in the wrapping balance becomes more 



important ||]. It is therefore desirable to have a better 
theoretical understanding of how physical parameters like 
bending stiffness, lateral tension, or adhesion strength 
control the shape of the complex and under which cir- 
cumstances complete wrapping ensues. This question has 
been partially addressed in a recent work that considers 
the adhesion of a cylindrical rod to a membrane and esti- 
mates the force required for unbinding j^. However, the 
description is limited to the regime of small wrapping, 
since it relies on a small gradient expansion of the elastic 
energy. In the present paper we determine the complete 
structural wrapping behavior of a spherical colloid by nu- 
merically solving the full nonlinear differential equations 
which describe the shape of an elastic fluid membrane 
under a prescribed lateral tension adhering to a spheri- 
cal particle with some given strength. These exact results 
are then compared to analytical calculations employing a 
small gradient approximation. We find that in the latter 
case the shape of the complex is only predicted correctly 
for very small degrees of wrapping. However, the phase 
boundary toward envelopment is still accurately repre- 
sented over a somewhat larger range of penetrations, es- 
sentially in the regime of low tension. 



II. GENERAL ENERGY CONSIDERATIONS 

We consider a spherical colloid of radius a adsorbed 
on a deformable surface, as depicted in Fig. [l[ Follow- 
ing Helfrich |l^ , the elastic energy associated with some 
membrane deformation is, per unit area, 



(1) 



with ci and C2 the local principal curvatures, k the bend- 
ing modulus, and a the lateral tension. From the two 
elastic constants we can construct a length, A, according 
to 



(2) 



For instance, with a tension of cr ~ 0.02dyn/cm (a value 
commonly found for cell membranes [pdl) and a typical 
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FIG. 1: Geometry of the wrapping complex and membrane 
parameterization. A membrane adheres partially and cylin- 
drically symmetrically to a spherical colloid of radius a with 
a degree of wrapping given by z = 1 — cos a. Due to the 
possibility of "overhangs" it is advantageous to parameterize 
the membrane by specifying the angle tp as a function of arc 
length s. The more direct choice of measuring the height h 
as a function of radial distance r remains however useful for 
the small gradient approximation. 



bending modulus of k ~ 20 (where /cbT is the ther- 
mal energy) we obtain A c± 64 nm (which, incidentally, 
is about the same size as many viral nucleocapsids [p^). 
Membrane deformations on a length scale smaller than 
A are mainly controlled by bending energy, while tension 
is predominant on larger scales. In this study, we will 
restrict ourselves to the bending-dominated low tension 
regime, a < A. We will see that even then a compara- 
tively low tension can play an important role. 

A natural parameter specifying the coverage of the col- 
loid by the membrane is the degree of wrapping z = 
1 — cos a. It ranges from z = 0, when the colloid just 
touches the surface, up to z = 2, in the fully enveloped 
state. The equilibrium value of z results from the adhe- 
sion of the colloid, driven by a contact energy per unit 
area, w, being balanced by the requirement to bend the 
membrane as well as the work of pulling excess membrane 
toward the wrapping site against the prescribed lateral 
tension a. 

Energies pertaining to the adhering part can be calcu- 
lated easily: The adhesion energy equals —2Tra'^zw, while 
bending and tension contributions are found to be Attzk 
and Tra^z^CT, respectively. The bending and tension en- 
ergies of the free part of the membrane are more difficult 
to evaluate, since they require the determination of the 
equilibrium membrane profile, a task that will be pursued 
below. Before solving the full variational problem, it is 
worth noting that the situation is greatly simplified for 
a membrane with no lateral tension. In this special case 
the equilibrium profile is quickly seen to be a catenoid, 
a minimal surface with zero mean curvature. Hence, the 
only energy contributions stem from the wrapped part 
of the membrane. We then find that, for (7 = 0, colloids 
do not adhere at low adhesion energy w < Wc — 2k/ a^, 
whereas full wrapping occurs above Wc, with no energy 



barrier to be overcome. 



III. NONLINEAR SHAPE EQUATIONS 

As a first approach, it would be tempting to approxi- 
mate the energy of the free section of the membrane by a 
phenomenological line energy. However, neither the rela- 
tion between the line tension constant and the membrane 
properties k and a would be known, nor is the implied 
dependency on the degree of wrapping correct. In order 
to draw indubitable conclusions, it is thus advisable to 
determine the exact membrane profile. 

The energy i^free of the free membrane is the surface 
integral over the local bending and tension contributions 
and is thus a functional of the shape. By solving the 
corresponding Euler-Lagrange equations one obtains the 
"ground state" profile and thereby its energy. Much work 
along these lines has for instance led to a detailed under- 
standing of vesicle conformations . In the present con- 
text, and using the parameterization indicated in Fig. |^, 
the shape equations can be written as 
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The dot indicates the derivative with respect to the arc 
length s along the membrane profile, and the p's are the 
"momenta" canonically conjugate to the coordinates ip, 
r, and h. These equations are integrated numerically for 
an asymptotically flat membrane that touches the colloid 
tangentially at a specific point of detachment. The de- 
termination of the appropriate boundary conditions far 
from the colloid raises a few subtle technical questions, 
and we refer the reader to Refs. ||l^ and jl^ for further 
details. From the numerical solution we finally obtain 
Efrcc ~ and thus the total energy - as a function of the 
degree of wrapping z. 

The structural phase diagram depicted in Fig. || follows 
directly. We find that the transition from the free to the 
partially wrapped state is continuous and occurs at the 
same value Wc = 2k jc? for any value of the tension. In 
contrast to that, an energy barrier separates the partially 
wrapped from the the fully enveloped state, rendering 
this transition discontinuous. The transition lines merge 
at a triple point (wc, 0), as expected. Our approach also 
allows us to determine the height of the energy barrier 
(see Fig. which surprisingly originates predominantly 
from tension. For cr = G.02dyn/cm (typical cellular ten- 
sion [0), a = 30 nm (capsid radius of Semliki Forest 
Virus, an often studied example) and k = 20 /cbT, the 
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FIG. 2: Structural wrapping phase diagram in the plane of 
adhesion constant w and lateral tension a. The bold solid line 
indicates the discontinuous transition between partially and 
fully wrapped, and the short dashed lines are the spinodals 
belonging to it. The fine dotted line o = w ~ Wc indicates 
where the fully enveloped state has zero energy. 
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FIG. 3: Height of the energy barrier iJbarricr as a function 
of reduced tension. The inset illustrates the shape of the 
function E[z) at the value of w where the transition from 
partial to full wrapping occurs and illustrates the definition 
of the concomitant energy barrier. 



barrier has the substantial value of £'bairicr ~ 22 k^T, 
showing that it cannot be overcome by thermal fluctua- 
tions. This point is also illustrated in Fig. |^ by showing 
the two spinodal lines at which the barrier vanishes. Inci- 
dentally, starting with a fully enveloped state and reduc- 
ing the adhesion w, the unwrapping-spinodal is crossed 
only once one is already in the unbound regime, i. e., the 
hysteresis would be so pronounced that one "skips" the 
partially wrapped region upon unbinding. 



The energy of the free section of the membrane van- 
ishes in the limit of full wrapping, z ^ 2, which is remi- 
niscent of the case of an ideal neck connecting two vesicles 
|l6| . The reason is essentially that as the neck contracts, 
the necessary membrane back-bending must occur on a 
length scale much smaller than A and is thus entirely 
dominated by the bending energy. But then the mem- 
brane will just assume the shape of a catenoid which, 
even though highly deformed, has zero mean curvature 
and thus does not cost any bending energy. 



IV. SMALL GRADIENT EXPANSION, 
RIGOROUS RESULTS 

The exact shape equations are nonlinear and can 
only be solved numerically. In order to get analytical 
information about the triple point at (wcO), we use a 
single-valued (Monge-) parameterization of the surface 
profile ft.(r), where r = {x,y) spans the flat reference 
plane at /i = 0. This representation does not allow for 
overhangs and therefore can only describe the first stages 
of the wrapping process. Assuming moreover that mem- 
brane deflections out of the horizontal remain small, i. e., 
I V/i| <C 1, the energy of the free part of the membrane is 



(4) 



This quadratic expansion is valid only for small deforma- 
tions but it has the advantage to render analytical cal- 
culations tractable. The equilibrium profile arises from 
the stationarity condition SEficc/dh = and satisfies the 
linear Euler-Lagrange equation 
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The general solution of eq. (||) is h{r) — hi + h2\n{r/\) + 
h3Ko{r/X) + hilojr/X), with Kq and Iq the modified 
Bessel functions Far from the colloid, the en- 

ergy density has to remain finite and integrable, so that 
h2 = hi = Q. The integration constants hi and h^ easily 
follow from the boundary conditions requiring the profile 
and its slope to be continuous at the point of detachment. 
The minimum energy of the functional in eq. (1) is then 
found to be 
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(6) 



where we introduced the abbreviation k — \J z{2 — z). 

The small gradient expansion and the exact membrane 
profiles are shown in Fig. ^ for different values of the de- 
tachment angle. As expected, the quality of the approx- 
imation declines as the degree of penetration increases 
(note that the parameterization itself fails much later, 
namely for z > 1). To further illustrate the discrepancy 
between the exact and the approximate solution, we plot 
m Fig. I the energy of the free part of the membrane. 
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The inset shows the energy difference between the small 
gradient expansion and the nonlinear result. We see that 
the approximation is fairly reasonable for degrees of pen- 
etration up to z ~ 0.1. Beyond this value, one has to 
solve the full shape equation in order to give an accurate 
description of the deformation energy. 

Nevertheless, the small-gradient expansion improves 
our understanding of the system, since it allows us to 
work out the asymptotic phase boundaries. We consider 
first the transition between the free and the partially 
wrapped state. For small degrees of penetration the total 
energy of the system can be expanded up to quadratic 
order in z, giving 
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where 7 — 0.5772 ... is the Euler-Mascheroni constant 
For sufficiently small z the quadratic term is negli- 
gible and wrapping sets in as soon as the prefactor of the 
linear term becomes negative. This confirms the numeri- 
cal finding that the initial wrapping transition is contin- 
uous and always takes place at w = Wc, no matter what 
the value of the tension is. 

Even though the final state of the second transition to- 
ward full envelopment cannot be described within a small 
gradient expansion, its energy is exactly known, since the 
free membrane does not contribute to it. The location 
of this second transition can therefore also be predicted 
by equating the energy of the partially wrapped state, 
eq. (|7|), with the energy of the fully enveloped state. 



E = 4^2 (cr 



Sttk. Since this yields a transcen- 



dental equation, the relation between a and w on the 
phase boundary cannot be expressed in terms of simple 
functions. However, the asymptotic behavior at zero ten- 
sion, which remains quite accurate over the whole regime 
of low tension, is surprisingly simple: We find that, up 
to a complicated but small logarithmic correction, the 
transition is given hy a = w — Wc (see Fig. 0). 



FIG. 5: Energy -Efroc of the free part of the membrane for the 
tension aa^ / k = 0.1 as a function of the degree of wrapping, 
z. Solid, dashed, and dotted curves correspond to the full 
nonlinear solution, the small gradient approximation, and its 
quadratic expansion, respectively. The inset shows the rela- 
tive error of the two small gradient expressions. 



Unfortunately, the concomitant energy barrier cannot 
be obtained in small gradient expansion. Even if the 
transition toward envelopment takes place at small z, 
where the approximation still works, the barrier occurs at 
large penetrations, z > 1 (see e. g. inset of Fig. ^). This 
is even beyond the range of applicability of the Monge 
parameterization h = h{x,y). 



V. DISCUSSION 

The combination of both numerical an analytical ap- 
proaches presented in this work has provided a complete 
description of the equilibrium wrapping behavior. The 
structure of the wrapping complex, the transition lines 
in the wrapping diagram, and the barriers of the discon- 
tinuous envelopment boundary can be predicted, and the 
low tension regime can even be treated analytically. One 
particular consequence of the wrapping scenario is that 
colloid engulfment is extremely sensitive to the particle 
size. As long as the tension is low enough {a < w/2 
will do), a large sphere is wrapped much more easily 
than a small one. This could be checked experimen- 
tally for instance by spreading a bilayer across a hole 
separating two compartments and adding a polydisperse 
colloidal solution to one side. If there exists some generic 
attraction between the colloids and the bilayer, a popula- 
tion of bilayer-coated colloids should emerge in the other 
compartment with a different polydispersity distribution: 
The frequency of large colloids is significantly enhanced, 
and too small ones do not occur at all. 

An important issue that we have not addressed here 
is the effect of thermally excited membrane fluctuations. 
At very low tension, the undulations might prevent the 
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beads from adhering to the membrane. An initial study 
of finite-size effects on the Helfrich repulsion shows that 
the (entropic) barrier to overcome in order to reach the 
surface increases with particle size ||l8| . This would pro- 
vide an interesting counterbalance to the tendency to en- 
gulf preferentially large spheres. 

The level of description in this work is that of generic 
physical mechanisms, but several examples mentioned in 
the introduction are taken from biology, where a host of 
different effects often occur at the same time. Notwith- 
standing this difficulty, we believe it to be worthwhile to 
analyze these examples - in particular, viral budding - 
in terms of the mechanisms discussed here. For instance, 
the adhesion strength of viral capsids to membranes can 
be modified by changing the biochemistry or simply the 
chemical potential of the binding proteins which in most 
cases are responsible for the attachment. Moreover, the 
membrane tension of cells increases or decreases during 
times of heightened endo- or exocytosis activity, respec- 



tively, and is generally under active control of the cell 
Iprf . Thus, processes like viral budding and its later 
counterpart of unwrapping (after infection of a differ- 
ent host cell and fusion with its membrane) invariably 
involve a movement in the structural phase diagram of 
Fig. H along both axes. The concomitant "phase changes" 
predicted that way thus offer an understanding of the ob- 
served (un)wrapping processes in terms of the physical 
forces driving them. 



Acknowledgments 

We would like to thank Bill Gelbart, Shelly Tzlil, and 
Avinoam Ben-Shaul for many stimulating discussions on 
the subject. MD also gratefully acknowledges financial 
support of the German Science Foundation under grant 
De775/1-1. 



[1] Lodish H., Berk A., Zipursky S. L., Matsudaira P., Balti- 
more D., and Darnell J., Molecular Cell Biology, Freeman 
& Company, New York (2000). 

[2] Garoff H., Hewson R., and Opstelten D.-J. E., Microbiol. 
Mol. Biol. Rev. 62, 1171 (1998). 

[3] Boussif O., Lezoualc'h P., Zanta M. A., Mergny M. D., 
Scherman D., Demeneix B., and Behr J. -P., Proc. Nat. 
Acad. Sci. 92, 7297 (1995). 

[4] Caspi A., Yeger O., Grosheva I., Bershadsky A. D., and 
Elbaum M., Biophys. J. 81, 1990 (2001). 

[5] Henderson R. M. and Oberleithner H., Am. J. Physiol. 
Renal Physiol. 278, F689 (2000). 

[6] Dietrich C, Angelova M. and Pouligny B., J. Phys. II 
(France) 7, 1651 (1997). 

[7] Koltover I., Radler J. O., and Safinya C. R., Phys. Rev. 
Lett. 82, 1991 (1999). 

[8] Deserno M. and Gelbart W. M., J. Phys. Chem. B 106, 



5543 (2002). 

[9] Boulbitch A., Europhys. Lett. 59, 910 (2002). 
[10] Helfrich W., Z. Naturforsch. 28c, 693 (1973). 
[11] Morris C. E. and Homann U., J. Membrane Biol. 179, 
79 (2001). 

[12] Baker T. S., Olson N. H., and Fuller S. D., Microbiol. 

Mol. Biol. Rev. 63, 862 (1999). 
[13] Seifert U., Adv. Phys. 46, 13 (1997). 
[14] Jiilicher F. and Seifert U., Phys. Rev. E, 49, 4728 (1994). 
[15] Deserno M., manuscript in preparation. 
[16] Fourcade B., Miao L., Rao M., Wortis M., and Zia R. K. 

P., Phys. Rev. E 49, 5276 (1994). 
[17] Abramowitz M. and Stegun I. A. (eds.). Handbook of 

mathematical functions, Dover, New York (1970). 
[18] Bickel T., submitted. 



